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Entanglement monotones, such as the concurrence, are useful tools to characterize quantum corre- 
lations in various physical systems. The computation of the concurrence involves, however, difficult 
optimizations and only for the simplest case of two qubits a closed formula was found by Wootters 
[Phys. Rev. Lett. 80, 2245 (1998)]. We show how this approach can be generalized, resulting in 
lower bounds on the concurrence for higher dimensional systems as well as for multipartite systems. 
We demonstrate that for certain families of states our results constitute the strongest bipartite en- 
tanglement criterion so far; moreover, they allow to recognize novel families of multiparticle bound 
entangled states. 

PACS numbers: 03.67.Mn, 03.65.Ud 



Introduction — Entanglement proved itself to be a fun- 
damental concept in physics, with applications spanning 
virtually all areas of quantum science: these include an- 
tipodal topics such as the black hole information paradox 
and industrial realizations of quantum cryptographic de- 
vices. By definition, entanglement between two or more 
particles is given by those quantum correlations, which 
cannot be created by local operations or classical commu- 
nication (LOCC). For the quantification of entanglement 
one can use so-called entanglement measures or entan- 
glement monotones - parameters that indeed are non- 
increasing under LOCC. The concurrence and the entan- 
glement of formation are important parameters of this 
kind Q. 

A central problem for the quantification of entangle- 
ment is the fact that nearly all entanglement monotones 
are extremely difficult to compute. Indeed, most defini- 
tions of entanglement monotones contain nontrivial op- 
timizations, such as the optimization over all possible 
LOCC protocols or the minimization over all possible 
decompositions of a given density matrix. This difficulty 
is an important issue for the application of monotones to 
real world problems or experiments. 

A milestone in the theory of entanglement measures 
was the derivation of a closed formula for the concurrence 
of two qubits by Wootters in 1998 @- In this work, it 
was shown how the minimization over all state decompo- 
sitions can be done for such a special case. Consequently, 
the Wootters' formula has lead to many applications of 
the concurrence, e.g. for characterizing phase transitions 
in spin models [3|. In the following years, the formula has 
been shown to work also for a special type of multipartite 
measures by Uhlmann Q • Furthermore, the concurrence 
could also be computed for some special states with high 
symmetries 

In the present paper, we generalize the idea of Woot- 



ters to compute lower bounds on the concurrence. Our 
methods work for higher dimensional bipartite systems 
as well as for multipartite systems. Compared with the 
large amount of research about lower bounds on entan- 
glement measures @, 0| our approach has substantial ad- 
vantages: for the bipartite case we discuss a family of 
bound entangled states and show that our result gives 
the strongest separability criterion so far; for the mul- 
tipartite examples, our estimates give the precise value 
of the multipartite concurrence, and allow to identify a 
novel and simple family of bound entangled states. Fi- 
nally, our approach can also be used to estimate other 
quantities besides the concurrence, which might be useful 
to deal with entanglement monotones based on antilinear 
operators and combs Q. It should be noted that lower 
bounds on the concurrence based on Wootters' idea have 
appeared in the literature before Q ; as we will see, how- 
ever, the existing approaches are fundamentally limited. 

Setting the stage — To start, let us recall the main def- 
initions. For a general m x n-dimensional bipartite pure 
quantum state qab = IV'XV'I on Ha ® Hb, the concur- 
rence [J Q can be defined as 

C(\i')) = ^(l-Tr^), (1) 

where qa = ^b(\^)JM>\) is the reduced density matrix 
of the first particle [10( . A pure state is separable if and 
only if its concurrence is zero. The above definitions are 
extended to mixed states via the so-called convex roof 
construction, 

C(e)= min r ft C(W), (2) 

where the minimization is meant as an optimization 
over all possible ensemble realizations g = '^2 i Pi\'4>i){'4>i\, 
where pi > and '^2 i Pi = 1. The decomposition at- 
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taining the minimum is said to be the optimal decompo- 
sition. Clearly, this is a difficult optimization problem, 
some estimates have been obtained with different tech- 
niques [1, 0]- 

The bipartite bound — For our approach, we first need 
to reformulate the definition of the concurrence. The 
pure state \ip) can be expressed in a product basis as 
|V') = Sj=i Sj=i V-'ij'Ki)- Furthermore, we can define on 
Ha the generators of the group SO{m) as L a = — 
There are m(m — l)/2 generators of this type, 
similarly, there are n{n — l)/2 generators Sp of SO(n) 
chiTLb- Then, a direct calculation for the ipij shows that 
one can express the concurrence as (see also Ref. 



C 2 (bA)) = 2(l-Tr^) 



a,/3 



\L a ®S p \r)\ 2 , (3) 



where \ip*) denotes the complex conjugation. In the fol- 
lowing, it is convenient to use a single index for the ma- 
trices L a (g) Sp and we define Jt = L a ® Sp, where the 
index t runs from 1 to N = [mn(m — l)(n — l)]/4. 

In order to formulate our bound, we first fix an integer 
k. We then choose a subset of indices t = {t±, ...,t k } C 
{1,...,N}, where we use the ordering t, < U+i. More- 
over, we can choose k complex numbers u = {u s } for 
which the absolute values are bounded via \u s \ < 1. Then, 
we consider the quantity 



A k (g, t, u) = max {0, A« - V A«„}; 



i>i 



(4) 



here the numbers X mn are the square roots of the eigen- 
values of 



X = g(^2u s J ts )g*(^2u* s J t! 



(5) 



8=1 



in non-increasing order. Alternatively, one can say that 

(i) 

the Xmn are the eigenvalues of the hcrmitcan matrix 



Before proving this theorem, let us discuss some of its 
implications. Eq. (J7J is a lower bound for the concurrence 
for any given choice of the u. In order to obtain a good 
bound, the set of the u has to be optimized for the given 
state g. Often this has to be done numerically, but we 
will also present examples, where a good choice of the Jt s 
is given analytically. 

Second, for the case of two qubits there is only one pos- 
sible generator, namely L a = Sp = |0)(1| — |1)(0| = ia y . 
This implies that the only possibility in Observation 1 
is k = N = 1, and then Eq. (JT)) reduces to the well 
known formula for the concurrence of mixed states. Of 
course, obtaining a closed formula for the concurrence is 
a significantly more advanced result as one has to prove 
in addition that equality holds. In Refs. 0, Q this has 
been achieved by writing down an explicit decomposi- 
tion. This is, however, beyond the scope of the present 
paper, we will only focus on the problem of deriving lower 
bounds. 

Finally, one should add that other researchers have ob- 
tained lower bounds on the concurrence by using the for- 
mulation of Eq. (J3j> and ideas similar to the original con- 
struction 0- In these works, the terms \(ip\L a <S> Sp\ip*}\ 2 
are estimated separately. A single observable L a ® Sp, 
however, acts on a 2 x 2 subspace only, and for these 
subspaces the criterion of the positivity of the partial 
transpose (PPT) is a necessary and sufficient criterion 
for entanglement This implies that the approaches 
in Refs. [7| can never detect weak forms of entanglement, 
such as bound entanglement which is not detected by the 
PPT criterion [12J . On the other side, Observation 1, rep- 



resents a very strong criterion for bound entanglement, 
as we will see below. 

Proof of Observation 1. First we prove that for a fixed 
k, and fixed vector t we have that 



y = JVe(E u Me*(E u MVe- (6) 
V s s 

For our given k, we consider the set of all possible t and 
choose for any of them a different vector u and compute 
the corresponding A k (g, t, u). This leads to (^f) numbers 
and for these we can state our first main result: 

Observation 1. Let g be a density matrix acting on 
an to x n-dimensional bipartite quantum system and con- 
sider for fixed k all the possible t and a possible choice 
of u as discussed above. Then, a lower bound on the 
concurrence is given by 



C{gf 



> 



N 



t 



A k (gXu)f 



(7) 



Especially, if g is separable then Ak(g,t,u) = for any 
choice of k, t and u. 



s=l 



UaJi.IV>*>!} >A k (g,t,u), (8) 



where the minimum is taken over all decompositions 

g = 52iPi\i > i)(i } i\- Lct ^ and \Xi) bc tn e eigenval- 
ues and the eigenvectors of g. It is known that any 
decomposition of g is connected to the eigenvalue de- 
composition via a unitary matrix {/y, namely one has 
y/P~iHi) = YlT=i U ij(y/^d\Xj)) 01- Therefore, we have 



^/PiPji^il Z)«=i u * J t s \ipj) = (UYU T )ij, where the ma- 

\fKXp{Xa\ E s =i u s J t 3 
Y T 



trix Y is defined by Y a p = ^/X a Xp(Xa\ l^ s =i u sJt B \Xp) 
Since the J k are symmetric, the matrix Y = Y T is com- 
plex and symmetric and we can use Takagi's factorization 
[14| to write Y = VDV T with a real diagonal matrix D. 
The entries of D are nonnegative and given by the square 
roots of the eigenvalues of YY^ . Then, following directly 
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the argumentation of Ref. Q we have: 

k 

= i E i[^^ T ]«i} > ^ - E A ™t 

i i>l 

= A k (g,t,u) (9) 



where A™ are the entries of D in decreasing order. These 
quantities are, however, nothing but the eigenvalues of X 
in Eq. Therefore, if a state g is separable then a de- 
composition into pure states without concurrence exists. 
Due to Eq. © all the mean values of J k vanish, which 
implies already that A k (g,t,u) = 0. 

It remains to show that Aj-(g,t,u) can give a lower 
bound on the concurrence also for entangled states. 
Suppose that g = Yli Pil^i) (^1 ^ s an optimal de- 
composition of g. Then C(g) = Y^iPiC{\ipi)) = 

J2iPi\/j2tLi \(A\Jt\^*)\ 2 ■ From tnc argumentation 
above, we know that for fixed k and t and fixed ti, t k 
the estimates A k (g,t,u) < \(tPi\u s J t ,\ipi)\ < 

EiPii:LiM\Jt.m\ hol d fe 

Finally, using the rule (J2j=i x j) 2 — kJ2j x ] an( l the 
Cauchy- Schwartz inequality we can directly estimate the 
right-hand side of Eq. ([7]) as: 



E[M^)l 2 <^(T)c 



(10) 



The details of this calculation are given in the Appendix 
Al. This concludes the proof of Observation 1. □ 
Before proceeding to the examples, let us discuss the 
properties of the concurrence that were used in the proof. 
The starting point was Eq. ([3]) and the only further re- 
quirement needed was that the fact that the Jt = Jj were 
symmetric 0. Moreover, if A t = -Af were antisym- 
metric, then one has for any state |('0|^4t|'f/>*)| 2 = 0, so 
restricting to symmetric J t can be done without loosing 
generality. In summary, the convex roof of any quantity 
E(\ip)), which can be written as 



E 2 (m = Y,±™ t \(mt\r) 



(ii) 



can be estimated with our methods: one can split each 
Mt in a symmetric and an antisymmetric part and esti- 
mate the contributions from the symmetric part. The 
fact that some of the coefficients mt can be negative 
does not matter: using the relation ^ t | (-i/ 7 1 I V'* ) 1 2 = 1 
(where the Gt form an orthonormal basis of the opera- 
tor space) one can rewrite E 2 (\ip)) as a sum with only 
positive coefficients minus a constant term. 

Bound entangled states as an example — In order to 
show that Observation 1 results in a stronger separability 
criterion than best methods that are currently known, we 
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FIG. 1. (Color online) Detecting entanglement in the 
Horodecki 3x3 bound entangled state with white noise. The 
criterion of Observation 1 (OBS1) is stronger than previously 
known criteria. For values of p smaller than the values given 
by SEP the states g a (p) are proven to be separable. See text 
for further details. 



consider the family of 3 x 3 bound entangled states intro- 
duced by P. Horodecki [l6[ . This family of states g^ H is 
not detected by PPT criterion, but is nevertheless entan- 
gled for any < a < 1 . The detailed form of these states 
is given in Appendix A2. We consider a mixture of these 
states with white noise, g a (p) = PQa H + (1 — p)l/9 an d 
ask for the minimal p, so that the entanglement in g a (p) 
is still detected. First, we use Observation 1 with the pur- 
pose of detecting entanglement and find the optimal Jt 
via numerical optimization. We finally compare our val- 
ues with the values obtained via different known criteria: 
the Zhang-Zhang-Zhang-Guo (ZZZG) criterion [l7|, the 
Ma and Bao (MB) criterion [18j, and the method based 
on symmetric extensions and semidefinite programming 
(SDP) [Tii[2(J. We also used also the algorithm proposed 
in Rcf. [21] to prove separability of quantum states. This 
allows to compute values of p, for which g a (p) is provably 
separable. 

The results are given in Fig. 1. One clearly sees that 
Observation 1 provides the best criterion, but the com- 
parison with the separability algorithm also suggests that 
Observation 1 does not detect all states. 

Estimating the multipartite concurrence — For simplic- 
ity, we only discuss the three particle case, but our results 
can be directly generalized to arbitrary TV-partite states. 
Let us consider a pure state in a d x d x d-system. 
Its concurrence is given by 



C T (m = 03-(Tre?+Tre! + Tre!)], 



(12) 



where g\ = Ti'23(p), etc. are the reduced one-particle 
states. For this definition, it directly follows that for pure 
states C r (|V>)) 2 = ^([C (1|23) (l^»] 2 + [C( 2 I 13 )(|^))] 2 + 
[C (3|12) (l^»] 2 ), where C^ 23 ) (|^)), etc. are the corre- 
sponding bipartite concurrences. This definition is ex- 
tended to mixed states via the convex roof construction. 
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Clearly, C T (g) = if and only if g is a fully separable 
state. 

A first possibility to estimate the multipartite concur- 
rence is to start with an estimate of the bipartite con- 
currence for each bipartition (as in Observation 1), and 
then estimate the total concurrence C T from it. This is 
indeed a viable way, in Appendix A3 we present and dis- 
cuss a corresponding theorem. The disadvantage of this 
approach is that there are states which are separable for 
any bipartition, but not fully separable |22j. For them, 
this method will not succeed, since all the bipartite con- 
currences vanish. 

To overcome this limitation, one should note that 
C T (\tp)) 2 is of the same structure as Eq. (fTT|) : we define 
the operators J t as before, but separately for any biparti- 
tion and write J^' 23 = L^t&Sp 3 and similarly for the other 
bipartitions. Then we have the expression C T (\ip)) 2 := 

IEji<</>k 1|2 V>i 2 + i<^ 2|1 V)i 2 + MJ? ll2 mn 

So we have to consider 
Af(^f) 



max(0,A« -$>« n ), (13) 



»>i 



where the Ami are the square roots of eigenvalues of 



v tot W 7 1|23 . T 2|13 . 7 3|12n * 

s=l 
k 

„ V"~V * 7 1 ! 23 i * T 2 |13 . * 7 3|12x /-..n 

X 2^^ J * S +v s J t! + w sJ t ; ) (!4) 

s=l 

in decreasing order. Here, x = (u, v, w) denotes a triple of 
complex vectors which are normalized as in Observation 
1 and t = {t\, ...,tk}- For this quantity we can state the 
following: 

Observation 2. For any arbitrary mixed state on "H® 
H&H and for every fixed k and for arbitrary x we have: 



N 



£[(a 



k\QXS)f] 



(15) 



A proof is given in the Appendix A4. 

Multipartite Examples — We will consider two simple 
examples for three qubits, but these already demonstrate 
two interesting points: first, they give an idea how the ob- 
servables Jt and the coefficients x can be chosen; second, 
it turns out that the entanglement criterion in Observa- 
tion 2 is very strong and allows to identify a novel family 
of bound entangled states. 

As the first example, we consider the three-qubit 
Greenberger-Horne-Zeilinger (GHZ) state \GHZ) = 
(| 000) + | lll))/v / 2 and mix it with white noise, g G {p) = 
p\GHZ)(GHZ\ + (1 - p)l/8. Then we take the single- 
qubit operator S (a) = |0)(1| - |1)(0| and the two-qubit 
operator = |00)(11| - |11){00| and from them we 

form the operators J l ^ k = (8 L^ k ^ for all three bi- 
partitions. Applying Observation 2 for the choice k = 1 



and U\ = V\ = w\ = 1, one finds already from a single 
term in the sum of Eq. (|15[) that the three-qubit concur- 
rence is bounded by 



(C T lg G ( P )]f>l(ll5 P -l]) 2 . 



(16) 



For p = 1, this reproduces exactly concurrence of the 
pure GHZ state. Moreover, this bound shows that the 
state g G {p) is entangled for p > 1/5. This means that 
Observation 2 provides a necessary and sufficient crite- 
rion for entanglement for the family of states g G (p), since 
it is known that for p < 1/5 these states are separa- 
ble [23|. In fact, Eq. (|16|) gives a linear lower bound 
on the convex function C T [g G (p)] and this bound coin- 
cides with the exact value on the points p = 1/5 and 
p = 1. This means that the bound equals the exact value 
on the whole interval p e [1/5; 1] and for them we have 
C T [g G (p)} = (l{5p-l})/V6. 

As the second example, we consider the three-qubit 
W state \W) = (|001) + [010) + |100))/\/3 mixed with 
white noise, g w (p) = p\W)(W\ + (l-p)t/8. In this case, 
we use again the operator = |0)(1| — |1)(0|, but for 
two qubits we use the = |00) (10| - 1 10) (00 1 and from 
them we form the operators ,P^ k = <3)L^ k \ Applying 
Observation 2 for k = 1 and u± = Vi = tui = 1, we find 
that (C T [g w (p)}) 2 > (l/96)[p(8 + V3)- Vt] 2 , especially, 
the state g w (p) is entangled for p > p s = \/3/ (8 + \/3) w 
0.178. 

This is a very remarkable value for several reasons. 
First, using the separability algorithm from Ref. [2l[, one 
can prove that the states g w (p) are fully separable for 
p < 0.177, giving strong evidence that Observation 2 
provides a necessary and sufficient criterion for the family 
of states g w (p). 

Second, these calculations show that the states g w (p) 
exhibit quite peculiar entanglement properties: one can 
directly check that for p < 3(8\/2 - 3)/119 « 0.2096 the 
states have a positive partial transpose for any biparti- 
tion, and using the separability algorithm plj one finds 
that for p < 0.2095 the states are indeed separable for 
any bipartition. Hence, for p <G [p s ; 0.2095] the states 
g w (p) are separable for any bipartition, but not fully 
separable. This implies that they are bound entangled: 
no entanglement can be distilled from them, even if two 
of the three parties join. It was known that such states 
exist [22[, however, the existing examples required a so- 
phisticated construction. It is surprising that the sim- 
ple family g w (p) includes bound entangled states and it 
underlines the power of our approach that these states 
can be detected with Observation 2. Finally, the bound 
entanglement in the family g w (p) can easily be gener- 
ated experimentally (contrary to other known examples 
of bound entangled states) since adding noise to a pure 
state is easy in practically any experimental implemen- 
tation. 

Conclusion — We have provided a general method to 
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bound entanglement monotones by extending in nontriv- 
ial way the original construction of Wootters , an ap- 
proach that works for both bipartite and multipartite 
concurrence. We leave open the problem of determining 
for which states our method gives the exact value of the 
concurrence. It would also be very interesting to broaden 
our approach to the general classification of invariants of 
quantum states 
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APPENDIX 



Al: Details of the calculation in the proof of Observation 1 



Let us write 1^ = \(tpi\Jt k IV'i )l- By making use of the rule (J2j=i x j) 2 — k^2j=i x j an d the Cauchy-Schwartz 
inequality we can estimate the following: 



£[A fc (^«)] 2 <£(l>£r 



-.(0 



s=l 
k 



E^ElE^J + 2 E^EE r l' ) E r ^ ) 

i<j t h — 1 m—1 



t s=l 
k 



<E? 2 E fc E( r i l) ) 2 + 2 E^ 



k 




\ t m=1 



^ & E^EE( r i 4) ) 2 + 2fc E^' 

i t s=l i<j 



X EE(4% EE( r 

\ t h=l \ t W=l 



k 2 /N\^ fi v, 2k 2 /N\ v-^ 



^©(E^ 



JV 



E(r 



\ *=i \ t=i 



A' 



Ec 1 ^) 5 



A 



\ E( r *) : 



□ 



A2: The family of bound entangled states 
The family of bound entangled states from Ref. fl6| are explicitly given by 
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These states have a non- negative partial transpose and are not distillable, but they are entangled for any < a < 1. 



A3: Estimates on the multipartite concurrence from bipartite bounds 

We define the operators as before, but separately for any bipartition. That is, we write J^' 23 = ® S^ 3 and sim- 
ilarly for the other bipartitions. Then, we obtain the expression C T (\ip)) 2 := \ Y.M^T^W)? + ^A^W)? + 
\{tp\jf 12 \ip*)\ 2 ]. Now, let us consider 

A^ |23 (£, t, u) = max {o, A« - ]T ^ A»„} , (18) 
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where the numbers A™ are the square roots of the eigenvalues of g( J2s=i "s^' 23 )^ ( S«=i "s^tt' 23 ) m non-increasing 
order and similarly for the other bipartitions. We can then state: 

Observation 3. Lef g be a density matrix on a tripartite dx dx d-system and consider for fixed k all the possible 
t and a possible choice of u as in Observation 1. Then, a lower bound on the multipartite concurrence is given by 

C» 2 > \^LT, {[Af\g,t,u)f + [Af u)] 2 + [A 3 J 12 (^u)] 2 }. (19) 

V k J £ 

Here, the coefficients u can be chosen separately for any t and any bipartition. Especially, if g is fully separable then 
A^ bc (g, t, u) = for any choice of k, t and u. 

Proof. First, one finds in the same way as in the bipartite case: 

^[A^ieXti)} 2 <[c 1123 (q)] 2 (20) 



and analogous bounds for the other bipartitions. It remains to bound C T (g) from the values of C^ 23 (g), C 2 \ 13 (g) and 
C 3 l 12 (f?). For that, let us assume that g = '^2 i Pi\'4 , i)('4 , i\ is an optimal decomposition when computing the convex 
roof of C T (g). Let us denote Q^^ bc ^ = \(ipi\J^ bc \'4 , i)\ 2 ■ Using the Cauchy-Schwartz inequality we have that 



[C^igtf + lC^igtf + lC^ig)} 
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* (£^/£eP |23) ) 2 + (£^/l>P |13) ) 2 + (E^/£eF |12) y 

i y k i y k i y k 

= £ {p 2 E(©i 4)(1|23) + ©i l)(2|13) + e« (3|12) } 



+2£W£ep |23 W^ 



i<j y k y k y k y k 



+^£ e i 4)(3ii2 y£^' )(3ii2) } 

< J2 {p 2 £(e< 4)(1|23) + ©i l)(2|13) + ei l)(3|12) } 



+ 2 £ j£ i i)(1|23) + i ?;)(2|13) + ©i ?;)(3|12) J£ ef 1|23) + ei' ;)(2|13) + e« (3|12) } 

•i<i y /c y fc 

= 2£>^(|^))] 2 = 2C^) 2 , (21) 



which proves the claim. □ 



A4: Proof of Observation 2 



First, as in Observation 1, we can prove that 



t ^ £.p<im Em? 3 + ^ 2 S 113 + ^ 3|12 m>i > Ar (e,r, 
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By denoting 6p |6c) = C |VI)| 2 , 
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This concludes the proof. 
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